Depletion of the Bose-Einstein condensate in Bose-Fermi mixtures 
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We describe the properties of a mixture of fermionic and bosonic atoms, as they are tuned across 
a Feshbach resonance associated with a fermionic molecular state. Provided the number of fermionic 
atoms exceeds the number of bosonic atoms, we argue that there is a critical detuning at which the 
Bose-Einstein condensate (BEC) is completely depleted. The phases on either side of this quantum 
phase transition can also be distinguished by the distinct Luttinger constraints on their Fermi 
surfaces. In both phases, the total volume enclosed by all Fermi surfaces is constrained by the total 
number of fermions. However, in the phase without the BEC, which has two Fermi surfaces, there 
is a second Luttinger constraint: the volume enclosed by one of the Fermi surfaces is constrained by 
the total number of bosons, so that the volumes enclosed by the two Fermi surfaces are separately 
conserved. The phase with the BEC may have one or two Fermi surfaces, but only their total volume 
is conserved. We obtain the phase diagram as a function of atomic parameters and temperature, 
and describe critical fluctuations in the vicinity of all transitions. We make quantitative predictions 
valid for the case of a narrow Feshbach resonance, but we expect the qualitative features we describe 
to be more generally applicable. As an aside, we point out intriguing connections between the BEC 
depletion transition and the transition to the fractionalized Fermi liquid in Kondo lattice models. 

PACS numbers: 03.75.Hh, 71.35.Lk, 51.30.+i, 64.60.-i 



I. INTRODUCTION 

The Feshbach resonance has emerged as a powerful tool 
in studying ultracold atoms in regimes of strong interac- 
tions. For two isolated atoms scattering off each other, 
the Feshbach resonance is a singularity in their scatter- 
ing length due to the coupling of the atomic states to a 
molecular bound statei^ The singularity (at v = 0) oc- 
curs as a function of the detuning v, which is a measure 
of the energy difference between the atomic and molec- 
ular states. The value of v can be varied by an applied 
magnetic field, and this effectively allows one to tune the 
strength of the atomic interactions. 

For systems in the thermodynamic limit, with a finite 
density of atoms, there is no singularity at v — 0. Nev- 
ertheless, the vicinity of v = is a regime of interesting 
many body effects. For a Feshbach resonance between 
two identical fermionic atoms, the many body ground 
state changes from a Bose-Einstein condensate (BEC) of 
molecules (v <C 0) to a Bardeen-Cooper-Schrieffer (BCS) 
superfluid descended from a Fermi gas of atoms (y 3> 0). 
It is important to note that there is no true fundamental 
distinction between the BEC and BCS states here, and 
so the two limits are connected by a smooth crossover. 
Recent experiments 3 ' 4-5 on 6 Li and 40 K atoms have suc- 
ceeded in observing the BEC of molecules. 

The consequences of the two-body Feshbach resonance 
are very different for other atomic statistics. For a Fesh- 
bach resonance between two identical bosonic atoms, it 
has been argued recently^ that there is indeed a sharp 
singularity, i. e. a quantum phase transition, in the many 
body system as a function of v. This singularity is not 
precisely at v = 0, but is shifted away from it; it is not 



directly a reflection of the singularity in the scattering 
length of two isolated atoms, but is a new many body ef- 
fect. Here, the two limiting states are a BEC of molecules 
(u -C 0) and a BEC of atoms (v >■ 0). Unlike the ferm- 
ionic case above, these two states cannot be connected 
smoothly to each other. The fundamental distinction be- 
tween these states becomes apparent upon examining the 
quantum numbers of the vortices in the condensate: the 
quantum of circulation differs by a factor of 2 in the two 
limits, being determined respectively by the mass of a 
molecule or of an atom. 

In the present paper, we will consider the remaining 
case of a mixture of two distinct types of atoms, one 
fermionic and the other bosonic. Mixtures of fermionic 
6 Li and bosonic 7 Li atoms were studied by Truscott et 
aim- and Schreck et al.^ and they succeeded in achiev- 
ing simultaneous quantum degeneracy in both species 
of atoms. Recently, Feshbach resonances have been ob- 
served between bosonic 23 Na and fermionic 6 Li atoms by 
Stan et al.^ and between bosonic 87 Rb and fermionic 
40 K atoms by Inouye et alM> So the time is clearly ap- 
propriate to examine the many body properties of such 
mixed Bose and Fermi gases across a Feshbach resonance, 
in which the fermionic molecule of the two unlike atoms 
can also reach quantum degeneracy. 

Our primary result is that such a mixture of fermionic 
and bosonic atoms also has a quantum phase transition. 
Again this transition is a many body effect, and does 
not occur precisely at v — 0. We will map out the phase 
diagram as functions of v, temperature (T), and the den- 
sities of the atoms (see Figs. ^ El El El and also describe 
the strong consequences of thermal and quantum fluctu- 
ations in the vicinity of the phase transition. It is also 
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possible that the mixture phase separates: this will be 
studied in the Appendix, where we determine the region 
of instability to phase separation (see Fig.EJ. 

As in the Fermi-Fermi and Bose-Bose cases, the exis- 
tence of a quantum phase transition for the Bose-Fcrmi 
case can be easily understood by characterizing the two 
limiting cases. For v 3> 0, there is a BEC of the bosonic 
atoms and a Fermi surface of the fermionic atoms. In 
contrast, for v <C 0, there is a Fermi surface of the ferm- 
ionic molecules. If the number of the fermionic (Nf) 
and bosonic (Nj,) atoms are unequal to each other, for 
v <C 0, there will also be some residual atoms which 
are not in molecules forming their own ground state: for 
Nf > Nf, the extra fermions will form a separate Fermi 
surface of atoms, while for Ni, > Nf, the extra bosons 
will form an atomic BEC. We note that for Nt > Nb, 
scanning the detuning v takes us between limits with and 
without an atomic BEC. Consequently there must be a 
critical detuning at which the atomic BEC is completely 
depleted, and all the bosonic atoms have been absorbed 
into molecules. 

A novel feature of this quantum phase transition is that 
it can be entirely characterized in terms of the Luttinger 
constraints on the Fermi surfaces. 

(i) Consider, first, the phase without the BEC with 
v <C 0. Here, there are 2 Fermi surfaces, one with 
Fermi surface excitations which are primarily the 
fermionic atoms, while the other has Fermi sur- 
face excitations which are primarily the fermionic 
molecules. We establish in Section [V] that this 
phase obeys two Luttinger theorems: the atomic 
Fermi surface encloses a volume associated with 
precisely Nf — Nf, states, while the molecular Fermi 
surface encloses precisely iVj states. 

(ii) Now consider the phase with the BEC. When the 
BEC is small, this phase retains 2 Fermi sur- 
faces, one primarily atomic and the other primarily 
molecular. However, now the volumes enclosed by 
these Fermi surfaces are not separately conserved; 
only the total volume enclosed by both Fermi sur- 
faces is required to contain Nf states. Eventually, 
for v 3> 0, the molecular Fermi surface disappears 
entirely, and only a single Fermi surface with Nf 
states remains. The disappearance of the molecu- 
lar Fermi surface (in the presence of a BEC) is a 
second quantum transition whose character we will 
also discuss briefly in Section IVT1 

This paper will determine the value of the critical v for 
the BEC depletion transition, and describe critical fluctu- 
ations in its vicinity. At T — 0, we will find in Section IVTl 
that this critical point is generically in the universality 
clas a 12 i 13 of the density-driven superfluid-insulator tran- 
sition with dynamic exponent z — 2. There is also an 
interesting quantum multi-critical point for Nf = 2Nb at 
which the BEC depletion quantum transition has a dif- 
ferent character: this we will also describe. At T > 0, 



the BEC depletion transition is in the universality class of 
the A-transition of He, and so will display similar critical 
singularities: a peak in the specific heat, and anomalies 
in transport coefficients. 

Our quantitative results are determined within a mean- 
field picture, whose applicability is restricted to the case 
of a 'narrow' Feshbach resonance, where the relevant cou- 
pling is sufficiently weak. We nonetheless expect our re- 
sults to be at least qualitatively applicable to the (exper- 
imentally more common) 'wide' Feshbach resonance. 

We also find an additional T = quantum phase tran- 
sition involving the disappearance of the molecular Fermi 
surface. As shown in Section IVTl this is described by a 
z = 2 critical theory of free fermions. 

While our work was in progress, we learnt of the work 
of Yabu et a/J-L who addressed some related issues, but 
only in the limit of infinitesimal coupling between the 
atomic and molecular degrees of freedom. We will note 
their limiting results in Section IlIII 

We now outline the contents of the body of the paper. 

First, in Section [Q] we define the model Hamiltonian 
that will be used throughout the rest of the paper. In 
Section ITTT1 we consider the limit of vanishing coupling, 
where a purely classical analysis can be used^ 

Section HVI finds the phase structure for finite coupling, 
treating quantum-mechanical effects using a mean-field 
approach. In Section we describe our results regard- 
ing Luttinger 's theorem for the system. In Section IVll 
the mean-field result of Section II VI is reproduced using 
a field-theoretical approach, which further allows us to 
characterize the critical properties of the transition. 

In Sections lVIII and lVIIII two corrections are calculated 
to the mean-field theory, which can be used to determine 
the validity of this approximation. In Section IVIII the 
two-loop corrections to the free energy are found, while in 
Section IVIIII higher orders in the coupling are included, 
within a low-density approximation. 

In the Appendix, we consider the stability of the sys- 
tem against separation into two regions with differing 
densities. It is shown that the system is indeed stable for 
a broad range of parameters. 

II. BASIC DEFINITIONS 

The system consists of bosonic atoms b and fermionic 
atoms / which combine to form fermionic molecules tp. 
The energy, relative to the chemical potential /i, is for 
the atoms 

^ = e /_ M / = _^_ /i / (2.1) 



and for the molecule 

=£-M* = i-M* + «', (2-3) 
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including the detuning v. The masses obey 



m and, because of the interaction, the chemical 



potentials are related by fi^ = ^ + /i b 
The grand Hamiltonian is 



n = 



(Zlflfk + e k blb k + ZMiPk) 
9 I (4 + k>fkh>+blftA+k') 

Jk.k' 



V> /,t„ 



-A 



k,k'.£ 



(2.4) 



where J k denotes J d 3 k/(2ir) 3 . 

We assume throughout that the fermion spin is po- 
larized along some direction, so that both / and 
are treated as spinless. The fourth term (in g) causes 
the bosonic and fermionic atoms to couple and form 
molecules, while the final term (in A) is an interaction 
between pairs of bosons. We omit the interaction be- 
tween fermions because the exclusion principle forbids 
s-wave scattering between identical fermions and we as- 
sume that the interaction between / and ip will be less 
important than the coupling g. 

Taking the dimensions of momentum and energy to be 
unity, dim[/c] = dim[_B] = 1, we have dim[?/>] = — ~ and 
the same for the operators b and /. (Throughout, we 
shall measure temperature, energy and frequency in the 
same units, so that Ti = fee = 1.) The coupling constants 
have dimensions dim[g] = — ^ and dim [A] = — 2. 

At temperature T = 1//3 ^ 0, we have six dimen- 
sionless parameters. First let Nf, be the total density of 
bosonic atoms, including those bound in molecules, and 
let Nf be the same for fermionic atoms. (We consider a 
unit volume, so that density is synonymous with num- 
ber.) In the absence of any fermions, the bosons would 
condense at a temperature 



To 



2tt 



a 



(2.5) 



We can take as dimensionless parameters T/To — (3o/(3, 
N f /N b , mf/m b , u/Tq, 7 2 /T and A 2 (m b ) 3 T , where 
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£1 

8n 



1 f b\ 3 / 2 

2m' m \ 



(2.6) 



In what follows, it will not usually be necessary to take 
account of the coupling between bosons given by the final 
term of (|2.4|l . Except within the condensed phase, which 
will be treated in Section II VI the only effect of A is a 
renormalization of the boson mass, which we assume has 
already been incorporated into the definition of m b . 



A. Physical units 

In order to relate these parameters to experimental val- 
ues, we may choose a unit of volume of 10~ 15 cm 3 , which 



gives the unit of momentum as roughly 10~ 27 kg- m/s. 
Taking the unit of mass to be 6 amu, corresponding to a 
lithium-6 atom, the unit of energy is roughly 7x 10 -10 eV 
or 8 fiK. 

For a Feshbach resonance, we assume the expression 1 ^ 
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27rab g ABA/i 



(2.7) 



where etb g is the background scattering length, AB is the 
width of the resonance and Afj, is the difference in mag- 
netic moments. Using the observed background scatter- 
ing length between lithium-6 and -7 of a^ g = 2.0 nm, 9 
we may estimate the coupling constant. Taking, for in- 
stance AB = 1 G, A/i = /ie, the Bohr magneton, we find 
g ~ 1 in our units. For a boson density Nj, = 10 15 cm" 3 
and mass mb = mf = 6 amu, the value g = 1 gives a 
dimensionless coupling of 7 2 /T"o — 5 x 10" 4 . 

While the width of the resonance used here, AB = 1 G, 
is sufficiently large that A/i AB > To, it is nonetheless 
somewhat smaller than typical experimental values. For 
our purposes, a more relevant measure of the resonance 
'width' is the lifetime of the molecule state in the vacuum 
(for v > 0). This is calculated in Section lYUIl where we 
show that it is determined by the constant 7. Since the 
relevant energies are on the order of To, the condition for 
a narrow resonance is that 7 2 /To 1. For the numeri- 
cal results throughout this paper, we will always remain 
in this narrow limit, which is analytically more accessi- 
ble. As noted above, we expect our results to be at least 
qualitatively applicable even for the wider Feshbach res- 
onances observed experimentally. 

Following Ref. 0, we take 



A 



2tt 



-abb 



(2.8) 



where for a^, the scattering length for the boson-boson 
interaction, we use abb = 0.27 nm, 9 giving 



A 2 (m b ) 3 T = 2 x 10" 



(2.9) 



The detuning v appearing in the molecular dispersion 
relation (|2.3|l is given byi^ 



v = A\l{B - Bo) 



(2.10) 



where Bo is the magnetic field at resonance and B is the 
applied field. 



III. THE LIMIT g 







The case of vanishing coupling, which can be addressed 
with a classical approach, has been considered by Yabu et 
alH (The results presented in this section produce Fig. 3 
of Ref. Q, which corresponds to our Fig. |3 below.) 

For simplicity, we restrict the analysis to zero temper- 
ature, but similar arguments can be made for nonzero 
temperatures. We call the two Fermi energies e an d 
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£q , and the corresponding wavenumbers k$ and fcg • At 
zero temperature, all bosons are at e b — and fermionic 
atoms or molecules must be added at their respective 
Fermi levels. 

The atomic Fermi surface (FS) vanishes when all the 
fermionic atoms are contained in molecules, so that 



** = (fti a JV» 1 /8 



(3.1) 



(The number of states within a unit sphere in momentum 
space is 1/67T 2 .) For this arrangement to be favorable en- 
ergetically, the molecular Fermi energy, eft, must remain 
below the lowest atomic energy level. The boundary of 
the phase without an atomic FS is therefore where 



1 



<N f ) 



2/3 



= 



(3.2) 



2m^ y " J/ (6tt 2 ) 2 / 3 
Similarly, the molecular FS vanishes at the point when 



1 



2mf 



(Nf) 



2/3 



(6tt 2 ) 2 / 3 



= 



(3.3) 



The atomic (molecular) FS is therefore only absent for 
negative (positive) detuning v. 

To find the boundary of the phase with a BEC, we 
must consider the depletion of the condensate. Bosons 
will take fermions and form molecules as long as their 
final energy is lower, i.e. < £q +0. The phase boundary 
is therefore where £q = £q , which gives 



2m 



j(N f -N b ) 



2/3 



2m^ 



2/3 



(6tt 2 ) 2 / 3 ' 



(3.4) 



where the wavenumbers have been determined from Nb 
and Nf, using the fact that there is no condensate. 

It should be noted that, in this limit, the coupling to 
fermionic atoms reduces the tendency of the bosons to 
condense. (The same is true at nonzero temperature.) 



IV. MEAN-FIELD THEORY 

It is possible to go beyond the classical analysis used for 
vanishing coupling, by using mean-field theory. We will 
present here two parallel developments, in this section 
and Section IVTI respectively. The first is based on single- 
particle quantum mechanics, using the mixing between 
the fermionic dispersion relations caused by the presence 
of a BEC. The second uses a field-theoretic approach and 
considers perturbative corrections to the bosonic propa- 
gator. The former has the advantage of giving a some- 
what clearer physical picture and leading more directly 
to thermodynamic results (such as the question of phase 
separation, considered in the Appendix), while the latter 
leads naturally to higher-order corrections. 

In the remainder of this section, we present the 
quantum-mechanical approach, starting from the Hamil- 
tonian l|2.4|l . First, in Section llV Al we make a mean-field 



approximation and diagonalize the new Hamiltonian. We 
then find the condition that a BEC should be energeti- 
cally favorable, within this approximation. 

Since the Hamiltonian is defined in the grand canonical 
ensemble, we must then relate the chemical potentials to 
the particle numbers, in Section flV Bl Within the mean- 
field approximation, it is sufficient to find this relation 
to order zero in the coupling, neglecting two-loop cor- 
rections to the free energy^ Later, in Section IVIII we 
determine the higher-order corrections. 

In Section llV CI we restrict our attention to the case of 
zero temperature, where transitions occur between states 
with different numbers of Fermi surfaces. We identify the 
positions of these transitions and present the full phase 
diagram for T = 0. 



A. Mean- field Hamiltonian 

Replacing the boson field in (|2.4|l by a real constant 
ip gives 



H 



d 3 k r 



ilftfk 



(2tt) 3 
- + Xip 4 , 
which can be diagonalized to 
d 3 k 



it V'lV'fc - g<p(4>tfk + fU>k) 

(4.1) 



(27T) 



(4.2) 



The dispersion relations for the mixed fermions F, are 



z 



)±±y(^-$) 2 +4 ff V, (4-3) 



with the choice that £^ > for all k. 

Since the mixing will cause the dispersion relations 
to separate, the total energy of the fermions is lowered 
by nonzero tp. This quantum-mechanical effect, in con- 
trast to the purely classical effect described in Section 
IIIII therefore acts to favor condensation. 

We must analyze the energetics to determine the point 
at which a condensate becomes favorable. The grand 
free energy $ is minimized at temperature 1/(3 by a 
Fermi-Dirac distribution of each fermionic species x (x € 
{F, W}). Ignoring the thermal distribution of bosons, 
which does not depend on tp, the total free energy is 



where 



+ \p A + 



E RX &)i ( 4 - 4 ) 

a;6{F,*} 



d 3 k 
(27T) 5 



In 1 



(4.5) 



The phase transition to a state with nonzero <p occurs 
when the quadratic coefficient changes sign, i.e. when 



1 rf 2 $ 



= 



(4.6) 



<p=0 
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Specifically, for negative A, nonzero ip is energetically 
favored, so that the condensed phase is stable. Using 
<|P|l . (Q)l and (Q|l . we find 



d 3 fc n F (rf)-n F (^) 



( 2 -) 3 g-tf 



(4.7) 



where rip is the Fermi-Dirac distribution function. The 
integral equation A — may be solved numerically. 



B. Particle numbers 

Since experiments are necessarily performed at fixed 
particle number, the expressions for the numbers in terms 
of the chemical potentials must be found. Particles of the 
species b, f and ip are not independently conserved, so 
the relevant quantities are Nf and Nf» the total numbers 
of fermionic and bosonic atoms, respectively (including 
those contained in molecules). 

As mentioned above, it is sufficient within mean-field 
theory to determine these numbers to order zero in the 
coupling. Since the species F and \& each contain one 
atomic fermion, we have 



N f 



(4.8) 



The number of bosons is (with ub the Bose-Einstein dis- 
tribution function) 

/rl 3 h 
7^M££) + n F ($)cosH k (4.9) 

+ n F (^)sin 2 e k ] , 

where the first term represents the condensate, the first 
term in the integrand is the thermal distribution of the 
bosons, and 9 k is the mixing angled 

When ip = 0, such as along the boundary to the phase 
without a BEC, the expression for the number of bosons 
simplifies to 



Nh 



d 3 k 



*B$)+n P (#) 



(4.10) 



To locate this phase boundary for fixed particle numbers, 
we must find the values of \rf and /i b which give the 
required numbers and also satisfy A = 0. (Of course, 
a third parameter must be tuned to its critical value to 
satisfy these three conditions simultaneously.) 

Results from such a procedure are displayed in Fig. ^ 
which shows the boundary for Nt/Nb = 1.11 as a func- 
tion of the detuning v and temperature T — 1/(3- The 
masses of the atoms are equal, — m b , and the solid 
line has dimensionless coupling 7 2 /T"o = 2.5 x 1CP 4 . For 
comparison, the case of vanishing coupling, treated in 
Section lTTTI is also shown with a dashed line. Both curves 
reach the value T — To, as in the case of free bosons, for 
v — > oo. when molecules cannot be formed. 




0.3 0.4 
Temperature T/T 

FIG. 1: Phase boundary with detuning v and temperature 
T, for fixed particle numbers Nf/Nt — 1.11. The dashed 
line has vanishing coupling and has been found with a purely 
classical analysis. The solid line has dimensionless coupling 
J 2 /To = 2.5 x 10~ 4 , and has been determined using the mean- 
field theory of Section ITVI For both, the condensed phase is 
on the left-hand side (for lower T) and labeled by (b) 7^ 0. 
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FIG. 2: Phase boundary with fermion number Nf and de- 
tuning v, for three different temperatures. The coupling is 
l 2 /To — 2.5 x 10 -4 and the masses are equal, m* = m b . The 
two phases are labeled as in Fig. with the condensed phase 
favored for higher detuning, lower fermion number and lower 
temperature. 



In Fig.|2the same phase boundary is shown on a graph 
of fermion number versus detuning, for three different 
temperature values. The solid line is at zero tempera- 
ture, T = 0, while the two dashed lines have nonzero 
temperatures. The coupling is 7 2 /Xb = 2.5 x 1CP 4 and 
the masses are equal, = m b . As expected, Bose con- 
densation is favored by lower temperatures, as in the case 
of an isolated Bose gas. 

It remains to be shown that the system is stable against 
separation into regions with different densities. It is 
shown in the Appendix that it is indeed stable for a large 
range of parameter values. 
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FIG. 3: The phase diagram at T = with dimensionless cou- 
plings (a) 7 2 /To = lO -6 , (b) 7 2 /7o = 2.5 x 1(T 4 and (c) 
/To = 2.0 x 10 -2 . The atomic masses have been taken to 
be equal, m* = m b and the coupling between bosons is given 
by X 2 (m b ) 3 T = 2 x 1CT 3 . The three distinct phases have, 
respectively, no Bose-Einstein condensate and two Fermi sur- 
faces (labeled '2 FS, no BEC), a condensate and two Fermi 
surfaces ('2 FS + BEC), and a condensate and a single Fermi 
surface ('1 FS + BEC). The dotted line indicates the fermion 
number at which Fig. 0] is plotted. 



C. Zero-temperature phases 

At T = 0, the Fermi-Dirac distribution function is re- 
placed by a unit step and all bosons occupy the lowest- 
energy state. As noted by Yabu et a/.flA the phase dia- 
gram can be further divided into a region with two Fermi 
surfaces and a region with a single Fermi surface. (We 
ignore the trivial case without any Fermi surfaces, which 
requires Nf — 0.) 

Except when the atomic numbers precisely match, 
Nf = Ni,, the case of a single surface can only occur 
when there is a BEC. In this case, tp, the expectation 
value of b, is given by the minimum of the free energy $ 
given in l|4.4|) . so that we must solve 



2//V + 4\ip 3 



E 

xef-F,*} 



dR x 
dtp 







(4.11) 



(excluding the root tp = 0). 

Following the choice that £jf > in l|4.3[) . the second 
Fermi surface disappears when £jfL — 0, making the 
Fermi wavenumber for F fermions vanish. For this to be 
the case, we require ^ > 0, jjft > v and 



gtp = \Jfj.f(n^ - v) 



(4.12) 



which should be solved simultaneously with (|4. 1 1|> . 

These expressions, along with the results in Section 
IIV Bl for the particle numbers, allow the complete zero- 
temperature phase diagram to be plotted. In Fig. the 
phase boundaries are shown on a graph of fermion num- 
ber against detuning, for equal atomic masses m* = m b . 
The three sets of boundaries have couplings (a) 7 2 /T"o = 



1(T 6 , (b) 7 2 /T = 2.5 x 10" 4 and (c) 7 2 /T = 2.0 x lCT 2 . 
(Since the dimensionless coupling depends on the fourth 
power of the coupling g appearing in the Hamiltonian, 
these large changes in j 2 /Tq in fact correspond to changes 
in g of only factors of 4 and 3 respectively. All of these 
coupling values are within the narrow resonance regime.) 
Throughout, we take X 2 (m b ) 3 T = 2 x 1(T 3 . 

The boundaries divide the diagram into three regions, 
depending on the presence of a condensate and the num- 
ber of Fermi surfaces. In the region labeled '2 FS, no 
BEC, the discriminant A is positive, so there is no BEC 
and two Fermi surfaces. In the region labeled '2 FS + 
BEC, A is negative and there is a condensate, as well 
as two Fermi surfaces. The lowermost region of the dia- 
gram, labeled '1 FS + BEC, has a condensate and only 
a single Fermi surface. 

In the limit of vanishing coupling (as in Ref. ^] and 
Section llllfl . the boundary between the regions with one 
and two Fermi surfaces extends down to the line Nf — 0. 
The region with a single Fermi surface is then divided 
into two, with the left-hand side having a Fermi surface 
of molecules and the right-hand side a Fermi surface of 
atoms. Including the quantum-mechanical effects, these 
two regions are no longer distinct, with the single Fermi 
surface crossing over from having a molecular character 
on one side (lower v) to having an atomic character on 
the other (higher v). 

This crossover is illustrated in Fig. where the ef- 
fective mass m* at the Fermi surface is plotted. The 
fermion number is set at Nf = O.lNb and the coupling 
is 7 2 /To = 2.5 x 10~ 4 , so that the system is within the 
phase with a single Fermi surface (of \& fermions). The 
effective mass is defined as 



in 



d 2 e k 



dk 2 



(4.13) 



For v <C 0, the Fermi surface has an essentially molecular 
character and m* ~ while for v 3> 0, it is atom-like, 



with 



In Fig. |SJ the Fermi wavenumbers of the two fcrmionic 
species are plotted, for coupling 7 2 /T = 2.5 x 10~ 4 and 
two different fermion numbers, Nf — §A*fc (solid lines) 
and Nf = hNb (dashed lines). In both the phase without 
a condensate (solid lines for v/T < 0.25) and the phase 
with a single Fermi surface (solid lines for v/T Q > 2.9, 
dashed lines for u/T < -0.65 and u/T > 1.3) the 
wavenumbers are constant, due to the fixed particle num- 
bers. Only in the phase with two Fermi surfaces and a 
BEC do the Fermi wavevectors change with detuning. 
(At the fermion number used in Fig. 0] the system stays 
in the phase with a single Fermi surface throughout and 
fe* = kf, feg = for all detunings.) 

We now turn our attention to the line dividing the 
phases 2 FS, no BEC and 1 FS + BEC in Fig. This 
boundary is horizontal and starts at the point where the 
three phases meet; in Section fV Al we will prove that this 
is at exactly Nb = Nf . At this transition, two changes oc- 
cur: both the second Fermi surface vanishes and the BEC 
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FIG. 4: The effective mass m* at the Fermi surface, with 
fermion number Nf — 0.1%, coupling j 2 /To = 2.5 x 10~ 4 
and equal atomic masses m' — m b . As can be seen from 
the dotted line in Fig. these parameters give a phase with 
a single Fermi surface. This surface changes from having a 
molecular character, with m* ~ 



character, 
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to having an atomic 
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FIG. 5: The Fermi wavenumbers for the two mixed species 
of fermions, ^ and F, with coupling 7 2 /To = 2.5 x 10~ 4 and 
equal atomic masses rrJ = m b . The solid lines have fermion 
number Nf = %Nb, while the dashed lines have Nf = \Nb. 
As can be seen in Fig. |3] the solid line goes between all three 
phases (at v /To — 0.25 and v /To — 2.9), while the dashed 
line goes from the phase with a single Fermi surface to that 
having two and back again (at v /To — —0.65 and v/Tq — 
1.3). The wavenumbers are measured in units of kf, the Fermi 
wavenumber for free fermions with number Nf. 



FIG. 6: The phase diagram with the fermion chemical po- 
tential [i s plotted on the vertical axis and the detuning v on 
the horizontal axis. Both have been scaled to ^ -v, which 
is held fixed throughout the plot. The boundary between 2 
FS, no BEC and 1 FS + BEC in Fig. expands into a new 
phase, labeled '1 FS, no BEC, within which there are only 
molecules, whose density is constant (both Nf and Nb remain 
fixed in the shaded region). The atomic masses are equal, 
mf — m b , and the couplings are j 2 /To — 2.5 x 10 -4 and 
A 2 (m 6 ) 3 r = 2 x 10 -3 . (Since the boson density is not fixed 
in this plot, the value of To appropriate to the phase 1 FS, no 
BEC has been used to define the dimensionless couplings.) 



tical axis. Throughout the plot fi^ 



v, and hence the 



Fermi wavenumber for the molecules, fcg , is held fixed. 
In the region where > 0, the essential features are un- 
changed, with the same three phases as shown in Fig. 
The boundary between the phases 2 FS, no BEC and 
1 FS + BEC, however, is seen to extend into an entire 
phase, labeled '1 FS, no BEC. In this region, there is no 
condensate and \J is negative, so that there is only one 
Fermi surface, of molecules. This entire phase therefore 
has Nf = Nb and collapses onto a single line in Fig. [3] 
Moreover, because is constant, Nt and Nb are both 
fixed in this phase. The situation in the shaded region of 
Fig. resembles that in the Mott insulator lobes in the 
phase diagram of the boson Hubbard model (see Ref. ^3 
and Chapter 10 of Ref. E3); at fixed ^ , the density of 
particles is insensitive to the variation in the chemical 
potential pJ . 



appears, as the line is crossed from above. Physically, this 
results from the fact that molecules are highly energet- 
ically favored in this region, so that as many molecules 
as possible are formed, and the residual atoms form their 
ground state. For Nf > Nb, these atoms are fermionic 
and form a Fermi surface, while for Nf < Nb, they are 
bosonic and form a condensate. Precisely at Nf = Nb, 
there are no residual atoms, so that there is no conden- 
sate and only a molecular Fermi surface. 

This situation is illustrated by Fig. El which shows 
the same phase diagram as Fig. |3 but with the chem- 
ical potential for the fermionic atoms, pJ , on the ver- 



V. LUTTINGER'S THEOREM 

All the ground states in our phase diagram in Fig. |3| 
contain Fermi surfaces. In Fig. El we presented the evo- 
lution of the Fermi wavevectors of these Fermi surfaces 
in our mean-field calculation. In the present section we 
will discuss general constraints that must be satisfied by 
these Fermi wavevectors which are valid to all orders in 
the interactions. (Throughout this section, we shall be 
concerned only with T = 0.) 

We will base our arguments upon the existence of the 
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Luttinger-Ward functional^ <I>lw, satisfying 
^ 5*lw[G'] 



SG' 



(5.1) 



G'=G 



where G" is a dummy variable, G is the actual full 
(thermal) Green function and £ is the full self-energy. 
(Throughout this section, we will mostly be concerned 
with the full Green functions, which we shall denote with 
the symbol G. When we make reference to the free Green 
function, this will be denoted Go.) 

Following Ref. ITsl it possible to construct the 
Luttinger-Ward functional non-perturbatively, starting 
from the partition function Z of the system. It can be 
shown straightforwardly that, treating the Green func- 
tion as a matrix in its momentum (and frequency) in- 
dices, any unitary transformation of the free Green func- 
tion, Go — * UGolJ^ 1 , that leaves Z invariant also leaves 
$lw invariant. 

A standard proof of Luttingcr's theoremiS for a system 
of interacting fermions makes use of the invariance of 
Z under a shift in the frequency appearing in the free 
propagator, u> — > u> + a. In our case $lw is a functional 
of the three Green functions, one for each species, and Z 
is invariant under a simultaneous shift in two of the three 
frequencies, i.e. 

$ hw [G^(ivi),Gf(iv 2 ),G b (iu>)} 
= $ LW [G^(wi), Gf(iv 2 - ia), G b (iu + ia)} 
= <S> VN [G^'{ivi+i{i),Gt{iv2),G b {u J j + iP)] (5.2) 



for any a and (3. 

To proceed further, it is useful to set ^ b 



' and 



consider derivatives of the grand energy with respect to 
\J and /i^ . The derivative with respect to ^ yields 



(/V) - (b f b) =N f -N t 



Each term on the left-hand side can be rewritten in terms 
of the full Green functions, giving 



N f -N b = - 



d 3 k dui 

W¥' 



G{{iu) + G\{iu>) . (5.4) 



(The change of sign of the / term results from the anti- 
commutation of fermion operators.) 

From now on the manipulations are standard^ We 
make use of the identity 



G(iu) = iGUu))— E(icj) - — hxGUu) , (5.5) 

OLO OUJ 



which results from the Dyson equation. The first equa- 
tion of H5.2fl gives, together with 15. 111 . 



d 3 kdui 

WW 



= . 
(5.6) 

Combining these two with (15.41) and integrating by parts 
gives 



f " 1 (2^) 4 du 



lnG{M + lnG£M 



(5.7) 



(with the boundary terms vanishing because G(' 
for \uj\ — > oo). 

The integral over to can be treated as a contour integra- 
tion and closed above, due to the factor e luJ ° . Changing 
the integration variable to z = iuj replaces this by an inte- 
gral surrounding the left half-plane. Since both of the full 
Green functions G b k (z) have all their non-analyticities 



b ■ 



(5.3) and zeroes on the line of real z, we can write this as 



N f - N b =i 



f d 3 k 


r° dz 


d ■ 


1 (2tt)3 . 


/-oo 2^ 


dz . 



In G[ (z + i0+) + In G\ (z + i0+) - In G f k (z + «T ) - In G b k (z + i0~ ) 



r 



(5.8) 



The integral of z can be performed trivially to give 
d 3 k 



N f -N b = i 



(2tt)< 



ln ^ ( '° +) I ln^ (i ° +) 



G r k (iO-) 



G b k (iO~) 



(5.9) 

Using the analyticity properties of the Green functions, 
this gives 

W 3 t r 

' " (-^ + s'/) + e(-^ + ^ 

(5.10) 

where is the unit step function and S' is the real part 
of the self-energy evaluated for to = 0. 



N f ~ N b = 



(2tt) 3 



First, we consider the phase with no BEC. Here there 
are necessarily two Fermi surfaces, and, as we will now 
show, the volumes of the two Fermi surfaces are sepa- 
rately constrained, independently of the interactions. 



By definition, the absence of a BEC requires that there 
be no bosonic quasiparticle excitations at or above the 
chemical potential, so that the second term in the brack- 
ets in (|5.10|l vanishes. (Note that this does not imply 
that (b^b) = 0, which is not the case beyond mean-field 
order.) This leaves the statement of Luttinger's theorem 
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for this case: 

The right-hand side of this expression is interpreted as 
the (reciprocal-space) volume of the atomic Fermi sur- 
face. 

A similar result follows by taking the derivative of the 
grand energy with respect to ^ , which gives 



(5.12) 



Going through the same manipulations as above leads to 



N h = 



d 3 k 
(2^)3 



e(-^ + s^)-e(-^ + Ej) 



(5.13) 

corresponding to (|5.10(l . Since we are in the phase with 
no BEC, this gives 



N b 



(5.14) 



We have therefore proved that there are two statements 
of Luttinger's theorem in the phase with two Fermi sur- 
faces and no BEC. One, l|5.11|l . states that the volume of 
the atomic Fermi surface is fixed by the difference in the 
numbers of atomic fermions and bosons, while the other, 
Q5.14[) . states that the volume of the molecular Fermi 
surface is fixed by the total number of bosonic atoms. 

Now let us extend our considerations to the phases 
with a BEC. In Fig. we observe that these phases may 
have either one or two Fermi surfaces. Here we show that 
Luttinger's theorem only demands that the total volume 
enclosed within both Fermi surfaces is determined by Nf, 
the volumes of the two Fermi surfaces (if present) are not 
constrained separately. 

In the presence of a BEC, it is no longer the case that 
the second term in the brackets vanishes in l|5.10[l and 
(|5.1iS|) . Instead, if we add these two equations, we arrive 
at 



Nt 



d A k 



r6( 



d A k 



©(-tf 



(2tt) 3 v sfc fe ' J (2tt) 3 

The two terms in this expression give the volumes of the 
two Fermi surfaces. We see that their sum is constrained 
to equal the number of fcrmionic atoms. 



the phases with and without condensates, since the con- 
densate vanishes as this line is approached from below. 
Secondly, the line that divides the regions with one and 
two Fermi surfaces is the point where the atomic Fermi 
surface vanishes. Where the two lines meet, we see both 
that H5.11fl is satisfied and that its right-hand side van- 
ishes. We therefore have Nb = Nf. 



VI. QUANTUM PHASE TRANSITIONS 

We now present an alternative analysis using the lan- 
guage of field theory. In Section lYl Al we reproduce the 
result that l|4.7|) determines the presence of the conden- 
sate. Then, in Section IVI Bl we determine the boson 
propagator near the BEC depletion transition. In Sec- 
tion [V^Sl we describe the critical field theories for the 
various transitions. 

The dimensionless Euclidean action corresponding to 
the Hamiltonian (|2.4|1 is 

g v q 

- p5 £ {^iU~v h v + bpf q -piJp) ■ (6.1) 



p.q 



The symbol p stands for k and uj n , and likewise q for £ 
and v m , where the Matsubara frequencies uj n (u m ) are 
even (odd). The summations over p (q) represent sums 
over u n (v m ) and integrals over the momentum k (£). We 
have also defined 



— V 



(Gp) 1 = -iu)„ 



•6 



(6.2) 



the inverse of the free Green function, G p , and similarly 
Sq. (In this section and the following, we will use the 
symbol G to denote the free Green function, contrary to 
the notation of Section IVI) 

We have omitted from the action the coupling term 
between pairs of bosons, since we will be interested in 
the region near the phase transition, where this term is 
not important. 

Integrating out both of the fermions and expanding the 
resulting coupling term to quadratic order in b and b, we 
find that the effective action for the bosons is 



S. 



(2)r 



1 



.« &b} = j}^b p E° p b p 

p 



^j2 G i G t+ P ~ b P b p ■ ( 6 - 3 ) 

p,q 



A. Multicritical point 

A simple application of our statements of Luttinger's 
theorem can be used to show that the multicritical point, 
where the three phases meet in Fig. |3 (and where four 
phases meet in Fig. [BJ, occurs at precisely Nt = Nf. 

Firstly, according to l|5.11|l . the volume of the atomic 
Fermi surface is given by Nf — Nb, as long as there is 
no BEC. This is therefore the case on the line dividing 



A. Mean-field approximation 

By replacing b by a real constant tp, we should arrive 
at the results of Section IIVI In this approximation, we 
have 
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so that the coefficient is 
A = V + .9 2 / -0^^G{(^m)Gf(iu m ). (6.5) 

The phase transition will occur when the coefficient A 
vanishes. 

We can represent (|6.5II by 



O 



(6.6) 



which appears as a self-energy diagram in the boson prop- 
agator, drawn as a dashed line. (The two solid lines rep- 
resent fermion propagators.) 

The Matsubara sum can be performed by replacing it 
by a contour integration, giving 



N _ _ 6 , 2 / dH n F tfl) - n F (#) 



(2-) 3 d - at 



(6.7) 



in agreement with l|4.7|l . 



B. Boson propagator 

By retaining the frequency dependence of the boson 
field, but again keeping only terms quadratic in b and 6, 
we can determine the form of the long-wavelength, low- 
frequency excitations. 

The effective boson propagator is, from (|6.3(l . the re- 
ciprocal of 



56 
—p 



q+p 



(6.8) 



-iuj n + i\ + g 2 



1 



dH 



P v -iv m + i{ -%{v m + uj n ) + at 



£+k 



which replaces 
So(iw n ) = -icj n - /i b + g 2 



For k — 0, this gives 

dH n F (#)-n F (#) 



where the result 



n F (a — iuj n ) — n F (a) , 



(6.9) 



(6.10) 



for u n a boson Matsubara frequency, has been used. For 
small uj n , we can expand to give 



A — iw r , 



1 1 



1-5 



dH n F (tf)-n F $) 



^ {4 -it) 
dH n F (Z[) - n F $) 



(6.11) 



(Note that, as required, the coefficient of lo 2 ^ is in fact 
positive.) 

The effective boson propagator (for k = 0) is then 



G b (kj n ) = 



(6.12) 



with 



Z = 



1-9' 



and 

IS(wn) = Z 



dH n F (g) - n F (Zf) 

^ (a! - at) 2 

dH n F (a[)-n F (at) 



(6.13) 



A - co 2 n g 2 



(27T) 



(al-atr 



(6.14) 

The integrals in the expressions for both Z and £q(w„) 
diverge at zero temperature if Nf = 2Nt so that the two 
Fermi wavenumbers coincide. For any other parameters 
the integrals are finite, and the effective propagator has 
the form (|6.12|) . As we discuss in the following subsec- 
tion, this distinction leads to different field theories for 
the BEC depletion quantum transition for these cases. 



C. Critical field theories 

First, we consider the BEC depletion quantum tran- 
sition. This is the transition between the 2 FS + BEC 
phase and the 2 FS, no BEC phase in Fig. |3J The same 
theory also applies to the transition between the 1 FS + 
BEC phase and the 1 FS, no BEC phase in Fig. El The 
low momentum modes of the b boson field clearly consti- 
tute an order parameter for this transition. The effective 
action for the renormalized b field near the critical point 
can be derived by integrating out the fermionic degrees 
of freedom, as already outlined in Section IVIBI 

For Nf 7^ 2Nb, the integration of the fermionic exci- 
tations is entirely free of infrared singularities: the dif- 
ferences in the two Fermi wavevectors implies that there 
are no low momentum fermionic particle-hole excitations 
at low energies. The resulting action for b contains only 
terms which are analytic in frequency and momentum, 
and has the following familiar form: 



S c [b] = J d 3 r J dr 



6 t|*_ _L bW 2 b + s \b\ i 



(6.15) 



Note that the b field has been rescaled by a factor \J Z 
from the b field in Section IVI Bl and that its mass has 
been replaced by the renormalized mass fh h . The action 
S c [b] describes a quantum phase transition with dynamic 
exponent z — 2 driven by tuning the coupling s. This 
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transition has been discussed previously in Ref. and 
in Chapter 11 of Ref. 0. The upper critical dimension 
is d = 2, above which the quartic coupling u is formally 
irrelevant. Nevertheless, the coupling u is important for 
the T > crossovers in the vicinity of the quantum crit- 
ical point: these are as presented in Ref. fl3l 

For Nf = 2Nb, there are low energy fermionic particle- 
hole excitations at zero momentum, and so the above 
procedure has to be reconsidered. Now there are non- 
analytic terms in the effective action for b, but these have 
a structure similar to that found by HertaStt for the onset 
of ferromagnetism in a Fermi liquid. Evaluating (|6.8|) for 
this case following Hertz, we now find the effective action 



S c [b] = I d 3 k I du;\b(k,uj)\ 2 
d 3 r I ,/r 



A- 2 



S |6| 2 + ||6| 4 



(6.16) 



The bare — itu term in the boson propagator is not in- 
cluded above because it is less relevant than the non- 
analytic term generated from the Fermi surface excita- 
tions. The critical properties of the z = 3 critical theory 
in Eq. (|6.16() have been described earlier by Hertz, and 
the T > crossovers by Millis2i (see also Chapter 12 of 
Ref. E| ■ 

Next, we consider the critical theory of the 2 FS + BEC 
to 1 FS + BEC transition in Fig. EJ The same theory also 
applies to the 2 FS, no BEC to 1 FS, no BEC transition 
in Fig. Here, a Fermi surface disappears as its Fermi 
wavevector vanishes. The critical theory is then the z = 2 
dilute Fermi gas theory already discussed in Chapter 11 
of Ref. El All interactions are irrelevant for the crit- 
ical properties, and the quantum-critical crossovers are 
merely those of a free Fermi gas. 

Finally, consider the multi-critical point, noted in Sec- 
tion IV Al where all phases in Fig. and Fig. meet. 
Here, both the b bosons and the / fermions are criti- 
cal. The critical theory is merely the direct sum of the 
z = 2 dilute Bose and Fermi theories mentioned above. 
All interactions between the critical / and b modes are 
formally irrelevant in three spatial dimensions. 



tion lVIIIl we will find a new expression for A by replacing 
the mean-field theory with a low-density approximation. 

In the remainder of this section, we will use S b from 
(|6.8() to determine the corrections to the grand free en- 
ergy $. We will show that these are negligible, pro- 
vided that the dimensionless coupling j 2 /Tq is suffi- 
ciently small. We therefore require a narrow Feshbach 
resonance for the results to be quantitatively accurate. 

By integrating the effective action over the boson field, 
we arrive at an expression for the partition function in- 
cluding Gaussian corrections, 



Z (2) 



(det5Q(det 
(det S b ) 

The grand free energy is then given by 



(7.1) 



$ = ^ + ^ + ^h|sJ + ^G^| , (7.2) 



using H6.8f) . where 



*§ = ±^E lnS « 



(7.3) 



is the grand free energy for the species x in the absence 
of coupling. [The plus (minus) sign applies to bosons 
(fermions) .] 

Taking a factor of S b out of the logarithm gives <&q, so 
that the correction to 4> is 



$ _ $ 



o 



This is the full expression for the Gaussian corrections; 
to estimate the size of these corrections, we will calculate 
the result to order g 2 . Dropping terms of higher order 
gives a correction to the free energy $ of 



A$ - % E G q G q+P G P 



(7.5) 



VII. GAUSSIAN CORRECTIONS 

In order to test the validity of the approximations 
made, we shall calculate two different corrections to the 
mean- field results of the preceding sections. First, in this 
section, we find the corrections to the grand free energy 
$ to order g 2 . These will result in corrections to the 
expressions found in Section IIVBI relating the chemical 
potentials to the particle numbers. Subsequently, in Sec- 



which can be represented diagrammatically as 
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Reinstating explicit momentum integrals and frequency sums, we have 

~ 2 r d 3 k d 3 £ 



G b K, k) G f (u m ,£) G^{v m +uj n ,£ + k) 



(2tt) 3 (2tt) 3 

Both Matsubara sums can be performed using contour integration, to give 



(7.7) 



A$ = g 2 



d 3 k r d 3 e 



n F (f/_ fc ) - n F (tf )1 U(^) - "B(e/ - tf_ fc ) 



(2tt) 3 7 (2tt) 3 



tb ,t/ _ t^ 

Sfc + Q-fc Q 



(7. 



r 



after a change of variables, £ — > £ — k. 



A. Renormalization of the detuning 

As it stands, the integral over k is in fact divergent. 
As |fc| tends to infinity (with \£\ finite), the second Bosc 
factor, n B (£f — %{_ k ), tends to —1. In the first bracket, 

np(^) remains finite, so the integrand tends to ~ l/k 2 
and the integral over k is linearly divergent. 

This divergence can be understood by considering the 
self-energy diagram 

b 




(7.9) 



which gives a correction to the detuning v linear in the 
cut-off momentum, 



vq - 9 



d 3 k 2m f m b 1 



(27T) 



jV> k 2 



(7.10) 



where vq is the 'bare' detuning that appears explicitly in 
the action. 

We use this expression to write v , which appears 
within in (|7.1|l . in terms of u, and then keep terms 
only of order g 2 . The renormalized expression for A<i> is 
then given by 



A$ = g 2 



d 3 k f d 3 l 



MtLk) - Mtf)] [»b($) - n B (£t - e/_ fc )] 2m/m b n F {$) 



(2tt) 3 J (2tt) 



tb ,t/ _ tV 



k 2 



(7.11) 



where the dispersion relation now involves the renormalized (physical) detuning v. (We have retained the same 
symbols for the new, renormalized quantities.) 

This expression can be simplified somewhat by performing a further change of variable, taking k — > k + (m b /m^)£, 
and also making use of the result 



We have finally 



A3- = 



d 3 k f d 3 l 



(2tt) 3 J {^) 3 _^ k 2_ v 

2m. J m° 



n F (x)n B (y -x) + n F (y)n B (x - y) = -n F (x)n F (y) 



n F (H f )n B (e) - n F (f> B (£ fe ) + n F (^)n F (^) - n F (^) 



(7.12) 



2m^m b v 
vr$ k 2 



(7-13) 



in which the energies and £ b should be evaluated Note that there is no singularity in the integral over k 



at the following momenta: 



in i|7.13[l . since the numerator also vanishes at the point 
where 



(7.14) 
(7.15) 

(7.16) (for v > 0). 



1*1 = 



2m/ m b 



(7.17) 
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7: Two-loop corrections to particle numbers, shown as 
a fraction of the total particle numbers to lowest order. The 
masses of the particles and their total numbers are the same 
as in Fig. Q Nf/Nb = 1.11. The dimensionless coupling is 
j 2 /To = 2.5 x 10 -4 . At each temperature, the corrections 
have been evaluated taking the detuning at its critical value, 
from Fig. 



The expression for A$ must be differentiated with re- 
spect to the chemical potentials to give the correction 
to the number of each species of particle. The resulting 
integral can then be performed numerically. 

The results of such a calculation are shown in Fig. [7| 
where we plot the corrections to the particle numbers. 
These have been divided by the total numbers evaluated 
using the results of Section llV Bl The atomic masses have 
been taken to be equal, as in Fig. ^ and, at each tem- 
perature value, the detuning takes on its critical value. 

We have taken for the dimensionless coupling 7 2 /T"o = 
2.5 x 1CP 4 , as in Fig.^ so that this corresponds to a nar- 
row Feshbach resonance. Since the correction everywhere 
is less than 5% of the total numbers, we expect that the 
mean-field results provide a good approximation for this 
case. The magnitude of the correction scales with g 2 cx 7, 
so that the quantitative predictions become less reliable 
for a broader resonance. 



whose amplitude is proportional to the density of atoms 
in the system. (A fermionic atom must be present in the 
system initially, in order to couple with the boson.) Di- 
agrams involving further loops with bosonic or fermionic 
atoms, such as 




(8.2) 



are proportional to higher powers of the density. 

There is, however, a correction to the molecular prop- 
agator of order zero in density coming from the process 
shown in 1(7.9(1 . which can take place in the vacuum. To 
make a consistent low-density approximation, we must 
therefore correct the molecular propagator using a Dyson 
equation 




where only the term in the bubble diagram 1)7. 9|) of order 
zero in density is to be included. 
The boson self-energy diagram 



1> 



o 



(8.4) 



should now be used to give a low-density approximation 
for the phase transition. 



VIII. LOW-DENSITY APPROXIMATION 



B. Calculations 



In this section, we develop a low-density approximation 
involving all orders in the coupling g. 



A. Diagrammatic description 

In the mean-field approximation, we have included in 
the boson self-energy such diagrams as 



Or 



(8.1) 



The Dyson equation (|8.3|) gives the relation between 
the reciprocals of the bare and full Green functions 



"<? "9 a 2-^i ^q-p^p 
p 



(8.5) 



compare ((6.811 . [The sign difference results from the 
fermion loop in 1)6.6(1 .] 

Including only corrections to order zero in density, 
by dropping the Bose- Einstein and Fermi-Dirac factors, 
gives 



d 3 k 



1 



2m^m b 
m^k 2 



(8.6) 
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where q = (£, v m ). (The second term in the parentheses 
comes from renormalizing the detuning v 1 as in Section 
IVII Al above.) The integral can be performed analytically, 
to give 

5* = + 2 7V /^ -v-iv m , (8.7) 

where 7, defined in l|2.6|) . has been used. 

This function can be continued to one that is analytic 
everywhere except along the real axis, by replacing iv m 
by z. In terms of z, the full Green function is 

Gf(z) = — . (8.8) 

z-g-2>yy/t;t-v-z 

Along the real axis, the square root has a branch cut for 
z > (f — v which corresponds to the continuum of free- 
atom excitations. For v < 0, G^ has a single pole at the 



real value 

z = $ - 2 7 2 (l - , (8.9) 

corresponding to the renormalized molecule. For v > 0, 
there are no poles, since the molecule has a finite lifetime, 
decaying into two atoms. 

C. Spectral representation 

These analytical properties are best summarized using 
the spectral representation for G^, 

Gf(z)= ^ dx^-. (8.10) 
The 'spectral density' is given by 



x - U + v 



2 V- *Z1 + 9(-,) V ; 7 " s(x-^-2^L/T^-l)) 



_(8T11 



where is the unit step function and 5 is the Dirac delta 
function. 



1. Weak-coupling limit 

It should be noted that, in the limit of small coupling, 
7^0, the first term of l|8.11|l involves the Lorentzian 
representation of the Dirac delta function, 



D. Effective action for bosons 

Using the modified spectral density for the molecule, 
(18.11(1 , we can compute the quadratic term in the effective 
action for the boson field. Following Section IVI Al the 
modified coefficient is 



A = -/i 



6 1 9 2 



/J3* 1 _ 



lim 



Tr(t 2 + e 2 ) 



S(t) 



.12) 



For 7 small enough, the first term of l(8.11|l has weight 
only near x = , where x — + v can be replaced by 
v. The limit of vanishing 7 is therefore given by 



1*3 



x - u + v 



so that the spectral density becomes, in this limit, 



.13) 



= S(x-Ct), (8-14) 

which is precisely the result for the bare molecule, used 
in Section IVll 



We can now use the spectral representation for both 
Green functions and then perform the Matsubara sum. 
Since the spectral density for the / atom has the form 



jUx)=8{x-d) , 



(8.16) 



the low-density approximation to the discriminant is 
given by 



d 3 e 

(2tt) 3 



• x G 



(8.17) 

with pf(x) given by (|8.11() . above. This integral can be 
performed numerically. 

Since A provides a correction to A, we must find some 
measure by which to determine the significance of this 
correction. A comparison with A is obviously not possi- 
ble, since this vanishes everywhere along the mean-field 
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curve in Fig.^ Instead, we shall find a lowest-order cor- 
rection to the critical detuning by evaluating 



8v 



A- A 



dA 
dv 



(8.18) 



0,N f ,N b 

(Finding the curve A = exactly would require a much 
I 



larger computational effort, since A takes considerably 
more time to evaluate than A. Instead, Sv provides the 
first step of a solution of A = by Newton's iterative 
method.) 

The partial derivative on the right-hand side of (|8.18() 
is given by 



dA 

dv 



N f ,N b 



dA 



dv 



dA 
d^J 



djj_ 

dv 



Nj,N b 



r 



dA 



A 1 



~dv~ 



(8.19) 
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FIG. 



8: The phase boundary, without (solid line) and with 
(dashed line) the correction to the detuning of Section fVIIII 
The solid line is the phase boundary as in Fig. using 
the same parameters (and dimensionless coupling 7 2 /Ib = 
2.5 x 10~ 4 ). The dashed line includes the correction Sv, from 
18. 181 . found using a low-density approximation. For clarity, 
this correction has been exaggerated by a factor of 3. Since 
the correction is negligible for these parameters, the mean- 
field results of Sections ll VI and IvTl are sufficient. 



where all the derivatives are to be taken at constant (3. 
To lowest order in the coupling, this can be replaced by 



dA 
dv 



N f ,N b 



|H +0 (<? 2 ) , (8.20) 

uv / Nf ,N b 



using 1)4. 7fl . The quantity on the right-hand side is the 
increase in the boson chemical potential needed to com- 
pensate an increase in the detuning, and keep the particle 
numbers unchanged. This can be evaluated using the ex- 
pressions given in Section llVBI 

The results of such a calculation are shown in Fig. [5] 
The solid line is the phase boundary shown in Fig. ^ 
using the same parameters. The dashed line is the same 
curve with the quantity 8v from 1)8. 18|) added. (For clar- 
ity, we have multiplied Sv by a factor of 3.) The small 
magnitude of the correction suggests that the results pre- 
sented in Sections llVl and IvTl are valid, for the parameters 
chosen. In this case, it is therefore unnecessary to use a 



low-density approximation; the mean-field result is suffi- 
cient. 



IX. CONCLUSIONS 

From a theoretical perspective, the main contribution 
of this paper is a description of a quantum phase tran- 
sition with an intimate connection to the Luttinger the- 
orem. On one side of the transition (see Fig. [3J, in the 
2 FS, no BEC phase, there are 2 Fermi surfaces with a 
separate Luttinger theorem for each Fermi surface. Re- 
markably, one Fermi surface is constrained by the total 
number of bosons Nb, while the other is controlled by 
Nf — Ni,, where Nt is total number of fermions. These 
two Luttinger theorems are consequences of two number 
operators that commute with the Hamiltonian: the op- 
erator f^f — tfb in Eq. I|5.3[) . and the operator f^f + tp'ip 
in Eq. i|5.12[) . On the other side of the quantum crit- 
ical point is the 2 FS + BEC phase. Here a Bose- 
Einstein condensate is present, and the condensate ef- 
fectively thwarts one of the Luttinger constraints. The 
single remaining Luttinger theorem demands only that 
the total volume enclosed by both Fermi surfaces is con- 
strained by Nf. We presented a theory for this transi- 
tion, along with phase diagrams as a function of system 
parameters. 

It is intriguing to note a connection between the above 
quantum phase transition and a seemingly disconnected, 
recent analysis of a quantum phase transition in a Kondo 
lattice model of the heavy fcrmion compounds^ This 
was a model of electrons occupying localized / orbitals 
(the fa electrons, where a is a spin index) interacting 
the itinerant electrons in the conduction band (the Co- 
electrons). A boson, b, was introduced to represent the 
hybridization between the orbitals. The connection be- 
tween the Kondo lattice model and the model of the 
present paper now becomes clear once we identify the Co- 
electrons with the molecular fcrmionic state ip, and the f a 
electrons with the / fermions. The Kondo lattice model 
also has two number constraints analogous to Eq. I|5.3|) 
and Eq. I|5.12|l . with one crucial difference: the first con- 
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straint is local rather than global, and applies separately 
on each lattice site. With this mapping, the heavy Fermi 
liquid FL state of Ref. ^ can be identified with the 2 FS 
+ BEC and the 1 FS + BEC phases. Further, the FL* 
phase of Ref. |22j is the analog of the present 2 FS, no 
BEC phase. The FL* phase also has two Luttinger theo- 
rems, one fixing the volume of the conduction band Fermi 
surface of electronic quasiparticles, and the other the vol- 
ume of the 'spinon' Fermi surface. The presence of a local 
rather than a global constraint implies that there is an 
additional gauge force that affects the spinon Fermi sur- 
face and the quantum critical fluctuations of the Kondo 
lattice. Such gauge forces are absent in our present con- 
siderations of Bose-Fermi mixtures, but, apart from this 
absence, there is a remarkable similarity to the FL-FL* 
transition in Kondo lattice models. 

On the experimental front, an obvious signature of 
the quantum phase transition in the Bose-Fermi mix- 
tures is in the evolution of the Bose-Einstein conden- 
sate. It would be interesting to scan the detuning and 
look for the disappearance of the condensate fraction at 
the lowest temperatures. The corresponding "superfluid- 
normal" transition should also survive at T > 0, where 
its signatures are similar to the A transition in He. 

A more dramatic, but experimentally less accessible, 
signature of the transition lies in the values of the Fermi 
wavevectors, as sketched in Fig. |SJ Measuring the Fermi 
wavevectors would allow detection of a Fermi surface con- 
strained by the number of bosons, and its eventual evolu- 
tion across the transition to a Fermi surface constrained 
by the total number of fermions. 

Finally, it should be noted that we have not addressed 
here the alternative of a paired state of fermions. Since 
we have dealt with spin-polarized fermions, s-wave pair- 
ing between the atoms is excluded, but p-wave pairing 
remains a possibility^^* There is also the more novel 
possibility for pairing between the fermionic atoms / 
and the molecules tp, which could be favorable when the 
two Fermi wavenumbers are approximately equal. This 
would then lead to condensation of a composite boson 
comprised of two fermionic atoms and one bosonic atom. 
We intend to investigate this possibility further in future 
work. 



APPENDIX: STABILITY AGAINST PHASE 
SEPARATION 

This section uses the mean- field results of Section IrvH 
The temperature will be taken as zero throughout. 



1. The compressibility matrix 

To establish the stability of the system against separa- 
tion into two coexisting fluids, we evaluate the compress- 
ibility matrix, defined by 



a/3 



<9 2 $ 
' dn a dijft ' 



(A.l) 



for a,{3(z {f,b}. 

We now define the (canonical) free energy F(Nf,Nb) 
by a Legendre transformation, 

F(N f ,N b ) = <S>(n f ,^)+[i f N f +[i b N b , (A.2) 

where Nf and Nf, are the total number of Fermi and 
Bose atoms, respectively. (Note that the full fermion 
and boson numbers, which are conserved by the Hamil- 
tonian, are used.) The compressibility matrix K' is then 
the inverse of the Hessian of F, so that complete stabil- 
ity against phase separation requires that K' be positive 
semidefinite. 

It is in fact easier to work with the matrix K a p, given 
by the same expression IIA.lt . but with a,f3 £ {f,ip}. 
This amounts to a simple (but not orthogonal) change of 
basis; it is sufficient (and necessary) for K 1 to be positive 
semidefinite that K be the same. 

We begin with l|4.4(l and 14.5|l and use (|4.11|) to de- 
termine the implicit dependence of (p on the chemical 
potentials. We must then take second derivatives with 
respect to the two chemical potentials to find the com- 
pressibility matrix. In the presence of a condensate, this 
leads to an expression 



K, 



a(3 



K 



(o) 

a/3 



A 



(A.3) 



where K^p is the matrix of second derivatives, evaluated 
at fixed ip and r a is a function whose form will not con- 
cern us here. 

The denominator of the second term is 
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y(j- + ^), (a.4) 



where 



and 



rri^mf dn/dk 



(A.5) 



(A.6) 
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When A goes through zero, the determinant of K di- 
verges, so that the Hessian of F becomes singular, signi- 
fying that one of its eigenvalues vanishes. This marks the 
onset of instability; we conclude that stability requires 
that A > 0. 

When there is no condensate, i.e. in the phase labeled 
'2 FS, no BEC in Fig. it is found that the system is 
always stable. 

2. Physical interpretation 

The obvious physical interpretation of A is that it 
represents the resultant interaction between the bosons, 
coming partly from the explicit term A in the Hamilto- 
nian 1)2. 4)1 and partly from the interaction induced by 
coupling to the fermions. This induced interaction can 
alternatively be found directly by continuing the expan- 
sion 1)6.3)1 to fourth order in b and b. 

A resultant interaction of the form I) A. 4)) is familiar 
from the case where the molecular degrees of freedom 
are not included explicitly in the Hamiltoniani2£i2ii This 
corresponds to our model for v ^> 0, when only virtual 
molecules are formed and the coupling term tjy fb in the 
Hamiltonian (|2.4|l can be replaced by a boson-fermion 
scattering of the form pb'bf. The induced interaction 
then comes from the diagram 




(A.7) 



which gives a term proportional to the density of states 
at the Fermi surface (at T = 0) . 

In this case, the induced interaction is always attrac- 
tive, as can be shown by a simple physical argument. 
For experimentally accessible parameters, however, it is 
not strong enough to overcome the intrinsic repulsion be- 
tween the bosons, so that the phase is stabled In our 
notation, the boson-fermion scattering is suppressed by 
a factor of so that the induced interaction falls off 
as \jv 2 . For v <C 0, a similar picture is obtained, with 
the atomic and molecular fermions exchanging roles. 

In the case of intermediate v, the induced interaction is 
no longer so heavily suppressed, but it is also no longer 
the case that it is always attractive. The physical pic- 
ture is clarified in this case by rewriting the action in 
1)6.1)1 in terms of the fermions F and '5 introduced in 
Section HV Al These fermions are defined so that there is 
no coupling term in the action linear in if = (6); instead, 
the lowest order interactions have the form F ip 2 F and 
Fip 4 F, and the same for '5. The former reproduces ex- 
actly the diagram (|A.7|) above, with / replaced by F and 
'J: physically this is a boson-fermion scattering induc- 
ing an attractive interaction between the bosons, as de- 
scribed above. This accounts for the final term in (|A.4(1 . 
Note that the exclusion principle requires the momenta of 




-2-10123 
Detuning v/Tq 



FIG. 9: The phase diagram at T = 0, as in Fig. 03 with 
couplings 7 2 /To = 2.0 x 10" 2 and A 2 (m 6 ) 3 T = 2 x 10" 3 . 
The other parameters, and the labels for the three phases, 
are the same as in Fig. |3] The region where the phase is 
unstable, as determined in the Appendix, is indicated. 

the two fermion lines to be exactly at the Fermi surface, 
leading to z x being evaluated at fcp. 
The term F ip A F produces the diagram 

F, 




which also represents an induced boson-boson interaction 
and accounts for the integral in (|A.4(I . Since Wk > 0, it is 
always repulsive and represents the fact that the fermion 
energy is lowered by a uniform distribution of bosons. 

3. Results 

The sign of the resultant interaction A must be cal- 
culated numerically to determine whether the system is 
indeed stable. Using the parameters from Fig. [31 stabil- 
ity is found everywhere within the plot for cases (a) and 
(b). In case (c), where the coupling g is larger relative 
to A, there is a region of the diagram where the phase is 
not stable; this is shown in Fig. EH 

For large \v\ the attractive coupling from the diagram 
(|A.7(I is suppressed by a factor 1/v 2 as described above, 
so that the system becomes stable. (The region for large 
negative v is not visible on this plot.) For intermediate 
values of \i/\, the induced coupling becomes larger than 
the intrinsic coupling, A, and it is the competition be- 
tween the two diagrams (|A.7|I and (|A.8|I that determines 
the stability. 

Stability is therefore favored by a higher Nf/Nt>, since 
this increases A;* and hence the phase space for the di- 
agram <|A.8|) . The other diagram, l)A.7)> . increases more 
slowly with fc* since the internal fermion lines are re- 
stricted to be at the Fermi surface. For intermediate \v\ 
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and very small Nf/Nb, on the order of 1CP 3 , the instrin- 
sic interaction once more dominates the induced and the 
system is stable. This region is too small to be seen in 
Fig.H 

An analysis similar to that carried out in Ref. |2fj could 
be performed to determine the stabilities of the alterna- 
tive, mixed phases. It should be noted, however, that, as 
can be seen in Fig. [5J the boundaries between the three 
phases are not disturbed at the parameters we have con- 
sidered. 

Furthermore, the analysis above shows that increasing 



the coupling g (or equivalently 7 2 /Tq) beyond the value 
used in Fig. [5] would increase the value of \u\ required 
for stability at small Nf/Nb (ie. extend the unstable re- 
gion to larger \i/\), but would not decrease the stability 
at intermediate \v\. This follows from the fact that the 
latter is determined by the competition between the two 
diagrams l|A.7|) and (|A.8|I . whose relative magnitude does 
not depend on g. We therefore expect that, for a broad 
Feshbach resonance, there remains a large region of sta- 
bility for intermediate values of \v\, similar to that in 
Fig. El 
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